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ABSTRACT 


This final report summarizes progresses made in the four proposed areas of study: 
(1) scattering model development for sparsely populated media, such as a forested 
area; (2) scattering model development for dense media, such as a sea ice medium or 
a snow covered terrain; (3) model development for randomly rough surfaces; and (4) 
design and conduct of basic scattering and attenuation experiments suitable for the 
verification of theoretical models. Progress during the course of this research has 
been detailed in semi-annual reports and forwarded to NASA. Hence, only the final 
results in each category are given in this report. 
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1. INTRODUCTION 


A scattering model has been developed for a deciduous and/or a coniferous forested 
area. The components of the forest consist of leaves, branches, trunks, and a rough 
ground surface. These components are modeled respectively as needles or elliptic 
discs with various ellipticities, dielectric cylinders of finite length, large vertical 
dielectric cylinders, and an irregular, finitely conducting ground surface. Typical 
predictions of this model are given in the next section and a detailed model 
development is given in Karam and Fung [1989]. Other related studies include 
permitting near field interaction between scatterers [Fung, Chen and Lee, 1986], leaf 
shape effects [ Karam et al, 1988], extinction coefficients in forested areas [Karam et 
al, 1988], scattering properties of vegetation components [Karam, Fung and Antar, 
1988], scattering from a layer of defoliated vegetation [ Karam and Fung, 1988], 
empirical modeling of a vegetation layer [Liu and Fung, 1988], and emissions from a 
multi-layered random medium [ Karam and Fung, 1986]. 

A new approach has been applied to develop a dense medium scattering model and 
applied to scattering from a snow covered surface [Fung and Chen, 1988], This 
approach is based on the Helmholtz integral representation for the scattered field in 
terms of the field inside the scatterer. A pair correlation function accounting for the 
spatial distribution of scatterers is used to describe the nature of the dense medium. 
The general behavior of this model as a function of volume fraction and the size of the 
scatterer is illustrated in Section 3. Its application to imaging is discussed in Fung, 
[1988]. 

A rough surface backscattering model valid over the entire frequency range has been 
developed for a randomly rough, perfectly conducting surface [Fung and Pan, 1987]. 
Its range of validity has been determined and shown to be valid for a moderately rough 
continuous surface [Pan and Fung, 1987]. The surface rms slope should be less than 
0.4. A summary of the predictions of this model is given in Section 4. 

Experimental studies have been carried out in the generation of known rough surfaces 
[Rochier et al, 1988] for scattering measurements, polarimetric microwave scattering 
measurements from known rough surfaces [Blanchard, Nance and Chen, 1988], 
attenuation and scattering measurements from a dense medium [Blanchard, Syed and 
Fung, 1988], and bistatic microwave measurements of simulated natural structures 
[Blanchard and Nance, 1988]. Results of these studies have been included in the 
semi-annual report for the period May 1, 1988 to October 31, 1988 and, hence, will not 
be repeated here. 
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2. Behaviors of a Scattering Model for a Forested Area 


In this forest scattering model four major scattering mechanisms are included: 
scattering by the ground surface attenuated by the vegetation layer; volume scattering 
from the crown layer including leaves and branches; trunk-surface interaction 
attenuated by the crown layer; and surface-crown layer interaction. In what follows 
illustrations of these scattering mechanisms are shown for incident plane wave with q 
polarization (q=Vj or hj) and in 0j direction (0j is measured from the normal to the 
canopy). 


2.1 Ground Surface Scattering 


i 



z=0 

Crown 



Figure 2.1 Scattering Geometry 

For the scattering geometry in Figure 2.1, the backscattering coefficient due to the 
surface scattering can be written as [Karam and Fung, 1989a] 

Opq(s) = L lp L 2 p Rp q L lq L 2q (2.1) 

Where p is the scattered field polarization (q=v s or h s ) and Rpq is the surface 
backscattering coefficient. Ljq and L 2 q are the crown and the trunk loss factors which 
have Hj and H 2 depth respectively 
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( 2 . 2 ) 


Li p = e -[kl P Hl sec e * 1 
L 2p = e" [k2 P H2Sec0il 

and k^p and k 2 p are the crown and trunk extinction coefficient respectively [Karam et 
al, 1988b; Karam and Fung, 1989c]. 

2.2 The Crown Scattering 

The crown backscattering coefficient can be written as 


a p q (c)=4tt cosQi (n L < | F Lpq 1 2 > + n b < | F bpq 1 2 >) (2.3) 

where n^ and n^ are the leaf and the branch number densities. F^pq and F^q are the 
scattering amplitude tensors for leaves and branches in the backward direction 
[Karam and Fung, 1988, 1989b]. The ensemble average < > is taken over the leaf and 
the branch orientations. 


2.3 The Crown-surface Scattering 

The crown-surface interaction backscattering coefficient is 

a pq (c <-»s)= a pq (c s) +c pq (s -» c) (2.4) 


The first term in (2.4) is due to scattering from the crown followed by scattering from 
the surface. The second term is due to scattering from the surface followed by 
scattering from the crown. The explicit contents of these terms are 



(c — ^ s) = 4 tc cos0j L lp L 2p 


^lp l-'lq 

^lq — ^lp . 


^ pp ^2p 


{n L < I F' Lpq | 2 > + n b < | F' bpq | 2 >} 


(2.5) 
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Gpq(s -^C)=4l COSGj | ^_ k L i lq [ L 2q R ',q L lq ^q 

f"L< I F'Lp, ! 2 >+ n b < I F'bpq I 2 >1 (2.6) 

where R'pp is the Fresnel reflectivity. F'Lpq and F'^pq are the scattering amplitudes 
in the specular direction [Karam and Fung, 1988, 1989; Karam et al, 1989b]. 

When p=q, (2.5) and (2.6) reduce to 

o" p (c -» s) = O 0 „(.s -»t). 4ltH, Lj p L| p 

tn L < I F Ln 1 2 > + n b < | Fpp, 1 2 >)R' pp (2.7) 

2.4 Trunk-surface Scattering 

The trunk-surface backscattering coefficient is equal to 

Op q (t <->s)= Gp q (t -> s) + ajq(s t) (2.8) 

The first term is due to scattering from the trunk followed by scattering from the 
surface. The second term is due to scattering from the surface followed by scattering 
from the trunk. Since the trunks have dimensions larger than the incident wavelength, 
the cross-polarized component vanishes and the like cross-polarized component can 
be written as 

a£p(ws) = aJ p (s-H) = 4 jc H 2 n t [L lp L 2p ] 2 | F' tpp 1 2 R' pp (2.9) 

where n t is the trunk number density. F'tpp * s t * ie trun ^ c scattering amplitude in the 
specular direction. 


2.5 Results and Discussions 

Figures 2-7 are presented here to illustrate the contribution of each term to the total 
backscattering coefficient and the relative values of the like and cross-polarized 
coefficients. In these figures, the branches are modeled by randomly oriented, finite- 
length cylinders, and the leaves are modeled by randomly oriented circular discs. The 
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cylinder and the disc orientations are described by the Eularian angles of orientation 
(a, P, y). Due to the scatterer symmetry, the third angle of orientation y has no effect 
on the scattering coefficient. 

The selected cylinder parameters are as follows: a (radius) = 2.5 cm, h (length/2) = 
20cm, 0° < a < 360°; 10° < P < 70°; Ej. = 12.4 - j4.9; Vf (volume fraction) = 0.1x10'^. 
The selected disc parameters are as follows: a (radius) = 2 cm, h (thickness/2) = 
0.125 mm; 0° < a < 360°; 60° < P < 90°; mg (gravimetric leaf moisture ) = 0.46. The 
trunks are modeled by vertical cylinders with parameters: a (radius) = 12cm; H 2 
(length) = 10 m; e r = 12.4 - j4.9; n t (number density per unit area) = 0.1 1 m'^. 

The figures illustrate the backscattering coefficients as a function of the incident 
angles at three different frequencies. 



Incidence Angle ( Deg ) 


+ total 
A soil 
e trunk-soil 
□ branch-soil 
■ branch 
x leaf-soil 
A leaf 


Figure 2.2 The VV backscattering coefficient components for a 
forested canopy at L band (F = 1.5 Ghz) 

Figures 2.2 and 2.3 show the backscattering coefficient at 1.5 GHz. The dominant 
contributing term at small angles of incidence is from the ground, and at larger angles 
of incidence it is the volume scattering due to branches which dominates. Other terms 
are less important with the choice of parameters indicated. Figure 2.3 shows that the 
cross-polarization is due mainly to branches at this frequency. 
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Backscattering Coefficient ( db ) Backscattering Coefficient (db) 



Figure 2.3 The like and cross backscattering coefficient from 
a forested canopy at L band (F = 1.5GHz) 



Figure 2.4 The VV backscattering coefficient components for a 
forested canopy at C band (F=6.5Ghz) 
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Figure 2.5 The like and cross backscattering coefficients for a 
forested canopy at L band (F=6.5 Ghz) 

Figures 2.4 and 2.5 show similar calculations for the backscattering coefficients at 6.5 
GHz. Here the dominant contribution comes from leaves, with contribution from 
branches as a secondary factor. When frequency is further increased to 10 GHz we 
see from Figures 2.6 and 2.7 that the dominant contribution is from leaves. 



Figure 2.6 The backscattering coefficient (vv) components for a 
forested canopy at X band (F=10 Ghz) 
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Incidence Angle ( deg) 

Figure 2.7 The like and cross backscattering coefficients for a 
forested canopy at X band (F= 10 ghz) 

3. Behaviors of a Dense Medium Scattering Model 

In this section the theoretical results of a model for dense medium scattering is 
summarized. To demonstrate the scattering properties of the model, the simple case 
of scattering from a half space medium with densely packed spherical scatterers is 
considered. 

3.1 The Backscattering Coefficient 

This dense medium scattering model is based on the Helmholtz integral formulation 
for the scattered field in terms of the field inside the scatterer. The pair correlation 
between scatterers used to characterize the dense medium effect comes in when 
ensemble averaging is performed on the scattered field product to obtain the average 
power. The backscattering coefficient is given as 

°pp= 4?t kj cos0 s f(u,v)l u _ k()Sin 0 5 v=0 (3.1) 

where 
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(3.2) 


f(u,v)= 


n 0 87i 2 l|3l 2 D 


lk' e + k' z l 2 (k ei cos0 t + k' ei ) 

p=i ^ ,,T 'wy 


8rc 


(3.3) 


In equation (3.3), vq is the volume of each particle; Tp is the p-polarized Fresnel 
transmission coefficient, f is the volume fraction of scatterers; Icq is the wave number in free 
space and y (^ice ” ^snow^(^ice ^^snow^* ^ equation (3.2), nQ is the number density 
of the particles, 0 t is the angle of transmission between the air and the medium with a plane 
interface; and k e is the effective wave number of the scattering medium. The definition of 
other variables in the expression of f(u,v) are, 

k’ e = V k e - u 2 - V 2 , k’ z =^k^- u 2 -v 2 , k’ ei = Im (k’J, k ei = Im{k e }. 

For the sake of simplicity, the effective wave number kg is evaluated in accordance with the 
empirical formula of snow derived by Hallikainen [1987]. Note that the solution given in 
(3.1) is for single incoherent scattering only. Multiple scattering can be included by the 
matrix doubling technique. The dense medium scattering phase matrix needed in the matrix 
doubling can also be developed in a similar manner. 

The quantity D in (3.2) which accounts for the dense medium effect is obtained from 
averaging the scattered power. The process of averaging generates the "pair correlation 
function p(^ x ,^y,^ z )" of the scatterers. 

D = 1 + n 0 Jp(^4 y ,y e- jk ^df (3.4) 


where 


* A A A 

k = x (k o sin0 + u) +y v -z (k er cos0 t + k' er ) 

^ A A A 

£ = X + y£y + Z £z- 

Note that k er and k’ er are the real part of k e and k' e respectively. 
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Assuming an exponential, isotropic correlation function p(r) = e' r ^, the dense medium 
factor D in the backscattered direction can be integrated as 


D = 1 + 


47tn 0 L 2 e " 2a/L f 


f 


IK 


1+(IK L) . 


2 IKa + 


V 


2 IKL 
1+ClK L) 2 ) 


cos(2 IK a) 


+ 

v 


2a | 1— (IK L) 2 ^ 
L 1+(IKL) 2 ; 


sin (2 IK a) 


(3.5) 


where "a" is the radius of the scatterers and IK=[4kg sin 2 0 + (k er cos0 t + k' er ) 2 ] 1/2 . 

3.2 Theoretical Results 

In this section the effect of volume fraction, particle size, frequency, liquid water 
content, and correlation length on the backscattering coefficient of snow are 
presented. 



Figure 3.1 The effect of volume fraction of scatterers on the 
backscattering coefficient of snow. 

In Figure 3.1, the hh polarized backscattering coefficient at 45 degrees incidence angle 
and 18 GHz is plotted versus increasing volume percent of the scatterers for two 
particle sizes. The correlation length is chosen to be three times the particle radius for 
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each case. It is seen that the effect of increasing volume percent is to increase the 
scattering for both particle sizes. In the modeling of scattering from snow using either 
the Rayleigh or the Mie scattering model under the assumption that the particles are 
far apart from each other, unreasonably large particle size must be used to fit the level 
of the data. It is anticipated that the effect of closely packed scatterers is to increase 
the scattering such that smaller particle size can be used to realize the data fit. The 
increasing trend of scattering coefficient versus volume fraction observed in Figure 3.1 
agrees with this effect. 

In Figure 3.2, the hh polarized backscattering coefficient versus frequency is plotted. 
The angle of incidence is 45 degrees, the volume fraction is 30%, the liquid water 
content is 0.5%. Two different particle radii, 0.25 and 0.5 mm, are plotted. The 
correlation length is again chosen to be three times the radius. It is seen that the 
effect of increasing frequency is to increase scattering. 
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Figure 3.2 The effect of frequency on the backscattering coefficient of snow. 

In Figure 3.3, the effect of liquid water content on the backscattering coefficient from 
snow is plotted for hh polarization at the angle of incidence of 45 degrees. The volume 
fraction of the scatterers is 30%, particle radius is 0.5 mm, and the correlation length is 
1.5 mm. The effect of increasing water content is to increase absorption. The effect of 
increasing absorption is to decrease the scattering albedo. Decreasing the albedo 
causes a decrease in the scattering coefficient. This phenomenon is clearly 
demonstrated in Figure 3.3 for two frequencies, 9 and 18 GHz. 
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Polarization=hh 



Figure 3.3 The effect of liquid water content on the 
backscattering coefficient of snow. 

In Figure 3.4, the effect of the correlation length of the pair correlation function on the 
backscattering coefficient is plotted. The angle of incidence is 45 degrees, volume 
fraction is 30%, liquid water content is 0.5%, and the particle radius is 0.5 mm. 
Scattering coefficients at four frequencies, 3.0, 7.5, 18.0 and 25.0 GHz, are plotted in 
Figure 3.4. 



Figure 3.4 The effect of the correlation length of pair correlation 
on the backscattering coefficient of snow. 
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It is seen that at low frequencies of 3.0 and 7.5 GHz, the effect of increasing 
correlation length is to increase the scattering coefficient. At 25 GHz the trend 
reverses. At 18 GHz the trend is to first increase and then decrease as the 
correlation length increases. Actually, the 3.0 and 7.5 GHz curves eventually 
decrease at a much larger value of the correlation length. Also, the 25 GHz curve 
increases at small values of the correlation length. The phenomenon of increasing at 
small correlation values and decreasing at larger values is common to all frequencies 
and can be explained as a resonance effect. Depending on the frequency, the 
resonance occurs at different values of the correlation length. 

4. A Summary of a Surface Scattering Model 

In this section the results of the backscattering model by Fung and Pan [1987] are 
summarized and applied to investigate the polarization effects in rough surface 
scattering. 


4.1 The Backscattering Coefficient 

The vertically and horizontally polarized backscattering coefficients denoted by 
and for a randomly rough perfectly conducting surface with standard deviation a 
and normalized roughness spectrum W(K) is 


aj v>hh =(2k 2 /cos 2 0)e _4k2a2cos20 y[(4k z a z cos z e) Il /n!]W w (2ksin0,O) 


2n\n, .-|vi/(n)/ 


±4 sin 2 0 e _3k2 ° 2cos20 Y [(2k 2 a 2 cos 2 0) m /rn!]W (m) (2ksin0,O) 


2 Q \in , , ivi A m )/ 


n=l 


+ 4sin 4 0e _2k2 ° 2cos20 > [(kVcos 2 0) m /rn!]W (m) (2ksin0,O) 


(4-1) 


where k is the wave number; 0 is the incidence angle; the minus sign in the middle 
term is for HH polarization; and W( n ) (U,V) is the roughness spectrum of the surface 

related to the n^ 1 power of the surface correlation function by the Fourier transform as 
follows: 


W <n) (U,V) = “ e- jU4 - j V (iC) d^dC, (n=l,2 ) 

Note that the first term in (4.1) is the standard Kirchhoff term. It should be the only 
term left when we take the high frequency limit. This, indeed, is the case since the 
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exponential decay factor in front of the second and the third terms are larger than the 
corresponding increases in the infinite sums as the wave number increases. More 
specifically, the sum in the second term cannot increase faster than expOk^o^ cos^0) 
and, similarly, the sum in the third term cannot increase faster than exp(2k2(j2 CO s^0) 
as the wave number increases. 

In the low frequency limit we need to take only the n = 1 term in (4.1) and approximate 
the exponential functions by unity yielding 


a? Vjh h = 8kV W(2ksin0,O) • \ 


2 2 

(l+sin 0) , for vv polarization 
cos 4 0 , for hh polarization 


(4.2) 


which is exactly the same as the results derived from the first-order small perturbation 
theory. 

As an example, for a Gaussian correlation function, 

p n (U)=exp[-(4 2 +; 2 )/L 2 ] 


(4.1) takes the form, 


aj v>hh = (2kL/cos0) 2 exp[-(2kacos0) 2 ] T [(kocos0) 2m /(m! rn)]e“ (kLsin0)2/m 


m=l 


{4 m *± 2 m sin 2 0 exp(kccos0) 2 + sin 4 0 exp[2(kacos0) 2 ] } 


(4.3) 


The cross-polarized backscattering coefficient denoted by a® h is [Fung and Pan, 1987] 


O j v =a« h =(8/itcos 2 0)e- 2k2 ° 2cos2e 




2_2 2 Q \in+n 

COS 0) 


m,n=l 


m! m 


J 


dudv 


[u 2 v 2 W (m) (ksin0+u,v) W (n) (ksin0-u,-v)] 

,222 
k -u -v 


(4.4) 


It is easy to see that (4.4) vanishes in the high frequency limit due to the exponential 
decay factor. In the low frequency region we set m=n=l and approximate exp(- 
2k^a^cos^0) by unity yielding 
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8 4 4 4 

G hv = cj vh = — k a cos 0 

K 


■J 


u 2 v 2 W(ksin0+u,v) W(ksin0— u,— v) 
(k 2 — u 2 — v 2 )cos 2 0 


dudv 


(4.5) 


This expression is in agreement with Rice [1951] and also Valenzuela [1967] if we let 
the relative dielectric constant approach infinity in Valenzuela [1967]. Note that the 
singularity, [cos^0 (k^-u^-v^)]"!, j n (4 4.) an d (4 5) j s d ue to assumption of a 
perfectly conducting surface. It comes from taking the limit as e r tends to infinity on 
the following factor 


(e r -l) (e r -sin 2 0) 


l z le r (k z -u z -v 2 ) 1/2 +(e r k z -u z -v z ) 1/i l 


[e r cos0+(e r -sin 2 0) 1/2 ] [cos0+(e r -sin 2 0) 1/2 ] 


2 2 2U/2.-2 


while keeping (k^-u^-v^)'^ greater than zero. The dielectric formulation given above 
[Valenzuela, 1967] does not have a singularity when (k^-u^-v^)"-^ ten( j s to ze ro. 
Thus, in actual evaluation we should use the dielectric expression with a large enough 
e r . Figure 4.1 illustrates the convergence of in e r . It shows that when the 

relative dielectric constant is larger than 10000, the surface acts like a perfectly 
conducting surface in its scattering behavior. 



Figure 4.1 The convergence of IEM with Cj. 
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4.2 Theoretical Results 


For ease of reference we shall refer to (4.1) and (4.4) as the Integral Equation Model 
(IEM) since it is based on the integral equation for the surface current. In Figures 4.2 
through 4.5 we show the dependence of the like and cross-polarization coefficients on 
the roughness parameter ka between 0.25 and 1.0 while kL is kept constant at 3.14. 
We purposely select an intermediate kL value since in the low and high frequency 
regions, scattering characteristics can be easily determined from the small 
perturbation and the Kirchhoff models. For the purpose of obtaining some reference, 
the predictions of the Kirchhoff and the first-order small perturbation models in like- 
polarization and the second-order perturbation in cross-polarization are also plotted on 
the same graphs. At ka = 0.25 IEM predictions agree well with the small 
perturbation results in cross-polarization but deviate at larger angles of incidence in 
like-polarizations because the kL value is larger than what is required by the 
perturbation model. The Kirchhoff predictions are correct only at normal incidence. 



kL=3.14 

ka=0.25 

IEM (HH) 

IEM (W) 

IEM (HV) 

■ SPM (HH) 

a spm (W) 

A SPM (HV) 

+ Kirchhoff 


Figure 4.2 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ka=0.25 and kL=3.14. 
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ka=0.50 

— IKM (HH) 

— - IEM (W) 

- IEM (HV) 

■ SP (HH) 

□ SP(W) 

A SP (HV) 

+ Kirchhoff 


Figure 4.3 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ka=0.50 and kL=3.14. 



kL=3.14 

kc=0.80 



IEM(HH) 
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IEM (W) 


IEM (HV) 

■ 

SP (HH) 

□ 

SP (W) 

A 

SP (HV) 

+ 

Kirchhoff 


Figure 4.4 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ko=0.80 and kL=3.14. 
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kL=3.14 

ka=LOO 

IEM (HH) 

IEM (W) 

IEM (HV) 

■ SP (HH) 

□ SP(W) 

A SP (HV) 

+ Kirchhoff 


70 

Figure 4.5 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ka=1.0 and kL=3.14. 

As ko increases to 0.5, IEM predictions differ significantly from the small perturbation 
results in like-polarizations since the latter is no longer valid. The cross-polarized 
returns have also begun to differ, but not significantly. This shows that the cross- 
polarized coefficient given by the second-order perturbation method has a wider range 
of validity in kL and ko values than the first-order like-polarized coefficients. As ko 
increases further to 0.8 and 1.0, the differences in cross-polarized scattering become 
more significant and it is clear that the second-order small perturbation model predicts 
a faster drop-off angular trend than the IEM. In all cases shown the Kirchhoff model 
predictions for like-polarizations continue to lie in between IEM predictions. The 
overall angular trends for both like and cross are to drop off slower as ko increases, 
and the level increase in the cross is much faster than the like. 

Figures 4.6 through 4.9 show the IEM model behavior when the roughness parameter 
kL is varied from 1 to 6 while ka is kept at 0.2. In these graphs it is seen that for ko = 
0.2 and kL = 1 IEM agrees well with the small perturbation model in both like and 
cross-polarizations. As kL increases, the separation between VV and HH 
polarizations decreases, angular trends drop off faster, and IEM approaches the 
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Kirchhoff model in like-polarization. Due to the choice of the small ko value, the 
difference between the cross-polarized coefficients is not very significant even when 
kL is 6. This again shows that the second-order perturbation model for the cross- 
polarization has a wider range of validity than the first-order like-polarized model. 
Theoretically, this is to be expected. It is especially interesting to note that the level 
of the cross-polarized returns rises as kL increases from 1 to 3, but it then decreases 
for further increases in kL to lower levels than when kL was 1. The initial increase in 
level is due to having a larger scale of roughness, while the decrease is due to having 
a significantly smaller slope. This decrease is consistent with the fact that a 
Kirchhoff- type surface has a low level of depolarization. 



> — i 
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kL=1.00 

kc=0.20 

— . IEM (HH) 

— _ IEM (W) 

............. IEM (HV) 

■ SPM(HH) 
□ SPM (W) 
A SPM (HV) 
+ Kirchhoff 


Incidence Angle 


Figure 4.6 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ko=0.2 and kL=1.0. 
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kL=3.00 

ka=0.20 



IEM (HH) 
IEM (W) 
IEM (HV) 
SPM (HH) 
SPM(W) 
SPM (HV) 
Kinchhoff 


Figure 4.7 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ka=0.2 and kL=3.0. 
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kL=4.50 

ka=0.20 

IEM (HH) 
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- IEM (HV) 
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SPM (HV) 
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Figure 4.8 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ko=0.2 and kL=4.5. 
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kL=6.00 

ka=0.20 

IEM (HH) 

- IEM (W) 

IEM (HV) 
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+ Kirchhoff 


Figure 4.9 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ka=0.2 and kL=6.0. 


In Figures 4.10 through 4.13 we repeat the same calculations as in Figures 4.6 to 4.9 
for ka=1.0. In this case differences in level and angular trends are evident between 
the Kirchhoff, the small perturbation, and the IEM in both like and cross-polarizations. 
As in Figures 4.6 to 4.9, an increase in kL leads to smaller surface slopes and 
curvatures thus permitting better agreement between IEM and the Kirchhoff model in 
like-polarization. While the second-order perturbation estimate of the cross-polarized 
return does differ in every case from IEM due to the choice of ko = 1, they are in better 
agreement when kL is larger because the surface is smoother. 
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Figure 4.10 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ko=1.0 and kL=1.0. 
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Figure 4.11 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are kc=1.0 and kL=3.0. 
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Figure 4.12 Comparison of the backscattering coefficient for both like and cross- 
polarizations among DEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are k0=l.O and kL=4.5. 
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Figure 4.13 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ko=1.0 and kL=6.0. 
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To show frequency dependence we illustrate in Figures 4.14 to 4.17 the like-polarized 
scattering properties by varying ko from 0.125 to 1 and kL from 0.75 to 6, a change in 
frequency by a factor of 8. We see a gradual decrease in the separation between the 
VV and HH polarizations, a fast increase in the level of the backscattering coefficients 
in the low frequency region, a gradually faster angular drop off, and the approach to the 
geometric optics solution at the high frequency end. In the meantime it is seen that 
both like and cross-polarized returns agree with the perturbation solution at the low 
frequency end. Then, as frequency increases, the like-polarized coefficients begin to 
deviate from the perturbation solution first, and then the cross-polarized coefficient. 
Note that while the polarized coefficients approach the Kirchhoff solution at the high 
frequency end just like the case when kL increases, the cross-polarized coefficient 
near vertical incidence appears to saturate at the high frequency end. The difference in 
the cross-polarized behavior between large kL and high frequency cases is due to the 
fact that the surface does not change in the study of the frequency response, while 
permitting kL to increase means that the surface slope is becoming smaller. 



kL=0.750 

ko=0.125 


x 

IEM (HH) 
IEM (W) 

**•••• 

IEM (HV) 

■ 

SPM (HH) 

□ 

SPM (W) 

▲ 

SPM (HV) 

+ 

Kirchchoff 


Figure 4.14 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are kc=0.125 and kL=0.75. 
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Figure 4.15 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff and small perturbation models. The normalized 
roughness parameters of the surface are ko=0.5 and kL=3.0. 
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Figure 4.16 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ko=0.75 and kL=4.50. 
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Figure 4.17 Comparison of the backscattering coefficient for both like and cross- 
polarizations among IEM, Kirchhoff, and small perturbation models. The normalized 
roughness parameters of the surface are ka=1.00 and kL=6.0. 


5. DISCUSSIONS 

Significant progresses have been made in the surface scattering modeling and a new 
approach in generating a known randomly rough surface has been developed. At this 
point in time sufficient data have not been gathered to validate the scattering model. 
This will be carried out in the future. 

In the past the Percus-Yevich’s pair correlation function, which was based on 
microscopic considerations, was used in the dense medium scattering model and, 
hence, the correlation function does not contain the correlation length parameter. An 
independent approach to develop a dense medium scattering model based on 
macroscopic considerations has been performed, which leads to a correlation function 
with its correlation length depending on scatterer dimensions. This model will be used 
to interpret field measurements. To verify the scattering it is necessary to perform 
experiments on well controlled dense media with a variable volume fraction. This is 
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yet to be done. The use of the field data alone does not provide enough cases and as 
detailed a ground truth to verify a scattering model. 

In conclusion, the most urgently needed results are complete sets of measurements 
from known surface and volume targets covering a range of target parameters over a 
wide range of frequencies, angles, and polarizations. We believe that scattering 
models must be compared with measurements from both data collected in the field, 
and data collected from known targets, to establish their applicability. 
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